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Abstract 


Tlie  problera  considered  here  is  that  of  the  transition  region 
separating  a  uniform  plasma  from  its  confining  magnetic  field.   The 
geometry  chosen  is  one -dimensional.   The  magnetic  field  profile 
within  the  boundary  layer  is  determined  by  the  current  distribution 
i.e.,  by  the  paths  of  the  ions  and  electrons.   The  paths  of  these 
particles  on  the  other  hand  are  determined  by  the  fields  in  which 
they  move.   An  exact,  self-consistent  solution  of  Maxwell's  equations 
and  the  equations  of  motion  of  the  particles  is  obtained.   A  detailed 
analysis  is  made  for  a  hydrogen  plasma.  These  results  are  then  com- 
pared with  those  of  a  plasma  whose  positive  constituent  is  the  proton 
and  whose  negative  constituent  has  a  mass  which  approaches  zero  and 
a  speed  which  approaches  infinity  in  such  a  way  that  the  energy  re- 
mains finite.   The  transition  layer  thickness  for  these  two  plasmas 
agree  to  one  part  in  10  .   It  is  shown  that  this  latter  case  is 
equivalent  to  treating  the  electrons  of  the  hydrogen  plasi-ia  in  the  guiding 
center  approximation. 


Introduction 

The  problem  of  the  transition  region  between  a  plasma  and  its 
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confining  magnetic  field  has  been  considered  by  a  nmriber  of  authors 

In  most  cases  the  emphasis  was  on  "single  orbit"  systems,  i.e.,  sys- 
tems in  which  the  orbit  of  all  ions  (electrons)  of  a  given  energy 
are  identical  except  for  a  simple  translation.   For  such  systems  an 
approximate  solution  is  possible  in  that  domain  where  the  ion  and 
electron  densities  are  nearly  equal.   An  exact  solution  can  be  ob- 
tained only  by  setting  the  electron  temperature  equal  to  zero  .   Our 
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An  exception  is  the  field  reversal  problem  considered  by  Harris. 
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olDjective  in  this  paper  is  to  consider  a  somewhat  different  physical 
system;  one  in  which  it  is  possible  (and  reasonable)  to  choose  the 
electron  density  identically  equal  to  the  ion  density.   For  this  system 
we  shall  be  able  to  exhibit  an  exact  solution  without  requiring  that 
the  electron  temperature  be  zero. 

This  paper  is  an  extension  of  an  earlier  work  by  the  author 

to  include  the  electron  motion.   The  problem  is  quite  similar  to 

8 
that  considered  by  Tonks  ;  the  essential  differences  being  that  we 

allow  the  particles  a  component  of  velocity  in  the  direction  of  the 

magnetic  field  and  that  the  solution  is  exact. 

Briefly  the  problem  may  be  described  as  follows.   All  quantities 
are  assumed  to  depend  only  on  a  single  variable,  x.   The  magnetic 
field  is  directed  everywhere  in  the  z-direction,  and  is  uniform  from 
X  =  -  00  to  X  =  0.   As  we  cross  the  plane  x  =  0  we  enter  into  the 
plasma  domain,  where  the  field  decreases  and  eventually  becomes  uniform 
again  at  x  =  +  o°.   The  paths  of  the  plasma  constituents  are  determined 
by  the  fields  (magnetic  and  electric)  in  which  they  move.   Tlie  fields 
are  in  turn  determined  by  the  paths  of  the  particles.   Our  problem 
then  is  to  solve^self  consistently.  Maxwell's  equations  and  the  equations 
of  motion  of  the  particles. 

These  conditions  are  not  sufficient  to  determine  a  unique  solution. 
Even  if  the  distribution  functions  should  be  specified  at  infinity  there 
is  still  considerable:  freedom  as  to  how  the  particles  are  put  into  the 
transition  region.   Tliis  lack  of  uniqueness  allows  us  to  impose  a 
number  of  additional  simplifying  restrictions  on  the  solution.   The 


most  important  of  these  is  that  of  neutrality.   Given  the  ion  dis- 
tribution function  (to  be  prescribed  presently)  we  require  that  the 
electron  distribution  function  be  such  that  the  electron  density  is 
everywhere  equal  to  the  ion  density.   This  condition  gives  the 
electron  distribution  function  as  the  solution  of  an  integral 
equation.   We  shall  also  require  that  all  particles  have  the  same 
energy  and  that  the  distribution  functions  depend  on  the  z-component 
of  velocity  only  through  the  energy.   With  these  restrictions  and  a 
suitable  choice  for  the  ion  distribution  function  the  solution  is 
uni  que . 

Analysis 
The  only  one  of  Maxwell's  equations  that  relates  to  this  problem 
is 

where  j  and  j  are  the  ion  and  electron  currents  respectively  and 
are  given  by 

j_^     =     1  e  /  y  f _!_  {x,x,y,i)   dx  dy  dz  (2) 

and  f  and  f  are  the  ion  and  electron  distribution  functions.   The 
equations  of  motion  may  be  integrated  once  to  give  three  constants 
of  the  motion.   These  are  the  energy  and  the  generalized  momenta  in 
the  y  and  z  directions 


-  h  - 


E   =  2  m ,  (x  +  y  +  z  j 


+ 


P 
+ 


=  m_^  y  +  e  A/c  (3 ) 


where  A(x)  is  the  vector  potential  defined  by 

B  =  dA/dx 

and  normalized  so  that 

A(0)  =  0  . 

Since  the  distribution  function  must  be  constant  along  the 
trajectories  of  the  particles  in  phase  space  (Liouville's  equation), 
the  equations  of  motion  of  the  particles  will  be  satisfied  by  choos- 
ing the  distribution  function  so  that  it  depends  only  on  the  constants 
of  motion  E^  P,  and  Q  .   Writing  f(E,P^Q)  for  f  (x,x,y,z)  when  ex- 
pressed as  a  function  of  E,P  and  Q  the  current  may  then  be  written 

j   =+2£   r(p+eA/c)   f   (E,P,Q)  ^fiy^^j  dE  dP  dQ       (k) 
±       -  '^+   J  ±  S(E,P,Q) 

where   S(x,y, z )/S(E^P^Q)   is  the  Jacobian  of  the  transformation  and  is 

'2 


a(E,P,Q)     -m3i-r±V2m^E-    (P  ".     f    f 


See  for  example  Ref.  3  and  k. 
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A  factor  of  two  has  been  inserted  in  (h)   to  take  account  of  the 
multiplicity  of  the  mapping  (x  -*  -  x).   It  is  clear  that  the  domain 
of  integration  must  be  limited  by  the  inequality 

m^x   =  2m_^E  -  (P  +^)2  .  q2  >  0  ,  (5) 

that  is^to  the  domain  inside  a  parabola  of  revolution  whose  axis  is 
parallel  to  the  E-axis  and  whose  apex  lies  at  the  point  (O^  +  ek/c,   0). 

Our  problem  may  now  be  stated  analytically.    Given  f  (E,P,Q)  , 
construct  an  f  satisfying  the  integral  equation  W,(x)  =  N  (x)  where 

P    r     f .  dE  dP  dQ 

^  =  4-  /      ,    -  ,  -  (6) 


+ 


and  is  the  ion  (electron)  number  density.   The  domain  of  integration 
is  again  given  by  (5).   The  currents  may  then  be  determined  from  (4) 
as  functions  of  A.   Substituting  these  currents  into  (l)  we  obtain 
A(x)  as  the  solution  of  an  ordinary  second  order  differential  equation. 
We  shall  take  for  our  distribution  function  f  the  following 

expression 

W  m^ 

S(E  -  E_^),     P  >  V2m^  E_^ 


U 


W  \/2m^E^ 


0,  P  <  J^~T~ 

\  +  + 

Such  a  choice  requires  l)  that  all  ions  have  the  same  speedy  2)  that 
the  distribution  be  isotropic  and  uniform  at  x  =  «>  where  the  density 
is  W  and  3)  that  there  be  no  ions  in  the  domain  x  .<  0  or  equlvalently 

00 

in  the  domain  A(x)  <  0  since  we  have  chosen  A(0)  =  0  and  we  assume  that 


A(x)  increases  monotonically  with  x  (as  may  be  verified  from  the 
solution).   This  last  condition  follows  immediately  from  (5)  and 
(6)  and  the  additional  simplifying  assumption  that  f  is  independent 
of  Q.   In  fact  we  may  show  that  for 

N_i_  =  0,         X  <  0 

it  is  both  necessary  and  sufficient  that 


f^  =  0  ,  P  <  \/  2m_^  E_^ 


To  show  this  we  note  that  since  all  ions  have  the  same  energy  we 
may  evaluate  E  in  (5)  at  E    Thus  the  domain  of  integration  for  the 


P-Q  variables  is  a  circle  of  radius  \/2iii     E  centered  at  (eA/c^O) 
(see  figure  (l));  a  cross  section  of  the  parabola  of  revolution 
mentioned  in  connection  with  (5).   Since  we  require  W  (x)  =  0  for 
all  X  <  0  and  since  f  is  nonnegative  it  must  be  zero  within  all  circles 
in  the  P-Q  plane  whose  centers  were  located  at  (eA/c^O)  with  eA/c  <  0. 
But  as  we  are  looking  for  a  distribution  function  which  is  independent 


of  Q  we  must  therefore  set  f  ,  =  0  for  all  P  <  V2m^  E_^.   Briefly  we 

may  say  that  if  f  =  0  in  the  half  plane  P  <  \/2m~E^  then  certainly 

f  =  0  within  the  circle  P^  +  Q^  =  2ni ,  E  ,.   In  the  first  case  f^ 
+  +  +  + 

depends  only  on  P  and  in  the  second  on  both  P  and  Q. 

¥e  now  set 

2 

=  g(P)  5  (E  -  E  )  (8) 


^rt/ 


2m  E 


-  6a 


♦  Q 


/2m+E^ 


P=  /2m  +  E  + 


Figure  1.   The  domain  of  integration  in  the  P-Q  plane  is  a  circle  centered 
at  (eA/c,0)  of  radius  \/2in  E    The  requirement  that  there  be  no  particles 
in  the  domain  x  <  0  is  satisfied  by  setting  the  distribution  function 


equal  to  zero  for  all  P<  \/2m  E  or  equivalently  by  using  the  shaded  area 
as  the  effective  domain  of  integration.  ' 
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a+r 

0 

a  <    0 

Woo 

r       /-^    - 

Nc„ 

N^ 

/        g(?)   dP 

J   a 

0  <  cc  <  2R 

+ 

~  2r 

u 

2R 

^ 

a  -r 

2R 

2R  <    a 

"Where  g(P)  is  to  be  determined  by  the  condition  of  charge  neutrality 
imder  the  restrictions  that  g(<»)  =  1  and  g(P)  =  0  for  P  <  >/2m  E 
so  that  there  may  be  no  electrons  in  the  domain  x  <  0.   For  con- 
venience we  have  set  P  =  -  P.   Equating  the  ion  and  electron  densities 
given  by  (6 )  we  have 


(9) 


where  we  have  set  a  =  eA/c;  r  =  y2m  E~  and  R  =  >y2m  E  To  solve 
this  equation  we  differentiate  both  sides  with  respect  to  a  and  ob- 
tain the  difference  equation 


(a+r)-g(a-r)  =  — 

R 


One  may  see  immediately  that  the  solution  is  given  by  g(P)  defined 
in  figure  (2).   At  any  point  P  =  a  the  difference  between  the 
function  r  units  to  the  right  and  r  units  to  the  left  of  a  is  every- 
where equal  to  r/R  unless  a   lies  outside  the  interval  0  to  2R  in 
which  case  it  is  zero.   We  also  see  that  for  such  a  solution  to  exit 
it  is  necessary  for  R  to  be  an  integral  multiple  of  r.   This  result 
is  not  unexpected.   It  states  in  effect  that  the  electrons  will  com- 
pletely shield  the  ions  if  the  electron  orbit  fits  precisely  into 
the  ion  orbit  without  any  overlap. 


0 

a  <  0 

1 

0  <  a  <  2  R 

0 

2  R  <  a 

Now  that  the  distribution  functions   f     and  f     are  known  we 
may  calculate  the  current  densities  from   (h) .      Tliey  are 


^o  1^    ^2   -    -    ) 


/a 


a 


E 


j      f^  -  n   +  1)    (^   +  n   +  1) 


'o     r 


(10) 


when  n  is  an  odd  integer  defined  by  the  inequalities 


nr-r<a<nr+r 


(11) 


In  effect  it  tells  which  step  the  vector  potential  a   is  on  (see 

figure  (2)).   Both  ion  and  electron  currents  are  zero  for  a   outside 

the  interval  0  to  2R,   The  constant  j   (  =  eW  R/^im  )  is  the  maximum 

O        00  '   +' 

ion  current.   We  see  that  if  the  electron  and  ion  energies  are  equal 
(and  from  here  on  we  shall  assume  this  to  be  the  case)  their  current 
will  also  be  equal  (see  figure  (3)).   If  these  currents  are  substituted 
into  (1)  and  both  sides  multiplied  by  da  we  may  integrate  once  to  ob- 
tain 


a 


i(da/d^)  -  i(da/dx) 


t 


kne 


J     (j+  +  j_)  da 


(12) 


^o   « 


\ocf   _  1  ^qf 
.  R    3  R  ■ 


-(n-l)r  +  (a-nr+  r)    (12) 
L  3 


5  ( (a-nr  r   +  r-^) 


8a  - 


tg(P) 


I 

4r/R 

3r/R 

2r/R 

r/R 


3r 


5r 


7r 


9r       2R 


Figure  2.      Definition  of  g(P)  for  R  =  5r.      The  definition  may  be  extended  in 
an  obvious  way  for  ajiy  multiple  of  r. 


2r        4r        6r         8r        2R 


Hgure  3.      Ion  and  electron  currents  for  R  =  5r. 
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Tliis  is  essentially  an  equation  of  pressure  balance.   Grad" 
has  shown  that 

(+) 


a 


1 


"XX 


j_^  (a)  da 


(13) 


0 


(+) 


"Where  p    is  the  indicated  component  of  the  ion  or  electron  stress 

XX 

tensor  and  defined  by  the  relation 


(+) 


P  =  ^  , 

XX  + 


f      (x,y,z)  d:-:  dy  dz 


=  4  r/ 


m 


2ni     E   -    (p   +  eA/c)^   -    Q^     f  ,(E,P,Q)   dE  dP  dQ 


Equation  (l3 )  follows  immediately  on  integration  by  parts  with  respect 
to  a.      We  see  then  that  (l2)  na.y  be  written  in  the  form 


B  /8jt  =  B/Srt  +  p    +  p 

O'  '  XX        XX 


which  states  that  the  total  pressure  is  a  constant  throughout  the  plasma 

and  equal  to  its  value  at  the  interface.   Mote  that  the  particle 

(+) 
pressure  at  the  interface,  p    (O)  is  zero. 

'     XX 


In  the  problem  considered  here  the  ion  and  electron  press"iires  are 
Jo 


+ 

3 
"XX 


R 


L  R    3   R 


{Ih) 


"xx     gj.2 


2/ sin-   l)r^  -f  r  (a  -  nr  +  r)  --( (a-nr  )^  +  r^ ) 


Here  again  n  is  defined  by  the  inequality  (8). 
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It  is  worthwhile  digressing  for  a  moment  to  note  that  all 
quantities  vary  only  within  the  range  0  <  a  <  2R,  and  are  constant 
outside  (see  (l),  (8)  and  (o)).   Prom  (ll)  it  follows  that  a 
(=  eA(x)/c)  is  finite  for  all  finite  x,  and  therefore  that  the 
transition  layer  is  finite.   Furthermore  one  may  give  the  following 
physical  interpretation  to  the  transition  layer  thicl-iness  in  terms 
of  the  particle  trajectories.   Returning  to  equations  (3 ),  which 
relate  the  particle  trajectory  to  the  constants  of  motion, it  is 
easily  seen  that  the  ion  which  just  reaches  the  plasma-vacuum 


interface  (P  =\/2m,  E  )  with  no  component  of  velocity  along  the 
magnetic  field  (Q  =  O)  will  range  from  a  =  0  at  the  left  face  of  the 


boundary  to  a  =  2\/2m  E  =  2R  at  the  right  face.   Thus  the  boundary 
layer  thickness  is  characterized  by  the  orbit  of  that  ion  with  maxi- 
mum radius  of  curvature  at  the  plasma-vacuum  interface  (see  figure  {h)). 
Returning  now  to  (12)  and  integrating  a  second  time  we  obtain  x 

as  a  function  of  a 

a.   r  , 

PIP  op  1  "^^ 

+  l/Rr      |_  -   (n-   l)r-^    +r      (a-nr   +r)--   ((a-nr)"'    +r^)    )l     da 


\  -2 


(15) 


where  ^  2 

"D  "D 

P     =     ^^     =     64«  j^  R/3  e  B^      -  (16) 

B 
o 

In    (l6)  we  have  used   (l2)   evaluated  at  a  =  2R. 


10a  - 


vacuum 


B  =B, 


m.   X 


Figure  k.     The  edges  of  the  transition  layer  "bound  the  trajectory  of  the  ion 
of  maximum  "diameter". 
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Equation  (15)  together  with  (l2)^ which  may  be  written  in  the 
form 

B(a)  =  B^  1^  +  i  P  ((a/R)^-  ^(a/R)^  +  Vrx-^ 

-(n-l)r^  +  r   (a-nr  +r)  -  -((a-nr)"^  +  r-^ ) 
^3  J 


(it; 


determine  the  magnetic  field  parametrically  as  a  function  of  x.  It 
must  he  remembered  in  (1T)  that  n  is  a  discontinuous  fimction  of  a. 

Although  the  integral  in  (15)  is  elliptic  in  form  it  is  in 
general  more  practical  to  do  the  integration  numerically.   There  is 
however  a  limiting  case  which  proves  to  be  of  particular  interest 
and  for  which  an  explicit  solution  is  possible.   We  shall  let  the 
electron  mass  go  to  zero  and  the  velocity  go  to  infinity  in  such  a 
way  that  the  energy  remains  constant  (and  equal  to  the  ion  energ;y). 
In  this  limit  we  see  that  r  approaches  zero  and  n  approaches  infinity 
in  such  a  way  that  the  product  nr  approaches  a   (see  the  inequality 
(11)).   The  electron  pressure  (equation  (l^))  becomes  simply 

p-   =  2/3  ^  nr  =  2/3  -  a  •  (18) 

XX      '   e  '    e 

The  integral  (15)  may  now  be  performed  explicitly  by  solving  a  cubic 
equation.   For  p  =  2"  the  solution  is 


^  =  ^o  1/^  V|3  _  51  ^  12 

X  =  ^  2(22)"^/'^  [f  (8l.80,cp)-  F(81.80,  66.27)] 


eB 
o 


i     =   x/22  tan^  9/2  -  3 
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where  F(k^9)  is  an  elliptic  integral  of  the  first  kind.   In  figxire 

(5)  we  have  plotted  the  magnetic  field  profile.   It  is  noteworthy 

that  the  slope  of  the  magnetic  field  is  not  zero  at  a  =  0  and  o;  =  2R 

despite  the  fact  that  j  and  j  are  zero  there  for  all  finite  r  at 

these  points  (see  (l)).   The  slope  therefore  approaches  its  limit 

discontinuously  as  r  approaches  zero. 

The  transition  layer  thickness  in  a  plasma  for  which  P  =  2, 

measured  in  units  of  the  diameter  of  an  ion  in  a  uniform  field  B  , 

o' 

is  1.17535.   There  are  two  averages  with  which  this  result  might  be 

compared.   The  first  is  the  average  diameter  of  an  ion  orbit  moving 

in  uniform  fields  B  and  B  and  the  second  the  diameter  of  an  ion 

o       00 

moving  in  a  vmiforra  field  which  is  the  average  of  B  and  B  .   These 
are  I.2O7  and  I.I72  respectively. 

For  comparison  we  also  carried  out  the  integration  (l2)  numeri- 
cally for  a  plasma  for  which  P  =  2'  and  R  =  i|-3r.   For  a  hydrogen 


plasma  R  =  /  —  r  =  ■!+2.8r.   It  is  not  possible  to  distinguish 
e 

graphically  between  this  case  and  the  previous  one  where  r  =  0.   In  fact 
the  transition  layer  thickness  here  is  I.I753O  and  differs  from  the 
former  result  only  in  the  sixth  significant  figure.   This  very  close 
agreement  becomes  somewhat  less  surprising  when  one  observes  that  if 
the  integral  (15)  is  expanded  in  a  power  series  in  r/R  the  term 
proportional  to  r/R  vanishes  and  the  difference  between  the  two  results 
must  be  at  most  of  order  (r/R)  . 


12a  - 
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It  should  be  recalled  that  the  slope  of  the  magnetic  field 
must  he  zero  at  both  ends  of  the  transition  layer  for  a  plasma 
with  finite  r/R.   That  this  is  the  case  for  our  simulated  hydrogen 
plasma  is  not  apparent  from  figure  (5).   It  was  not  possible  to 
indicate  on  the  graph  the  very  rapid  fluctuations  in  the  slope  of 
the  magnetic  field  (see  (l)  and  figure  (3)  with  R  =  43r).   The 
slope  varies  from  zero  to  a  maximum  and  back  to  zero  again  in  a 
distance  which  is  roughly  l/i+3'rd  of  the  transition  layer  thick- 
ness.  The  average  slope  over  this  range  is  that  indicated  in  the 
figure. 

Guiding  Center  Approximation 

We  have  seen  that  the  agreement  between  the  "exact"  results 
and  the  results  obtained  by  setting  the  electron  mass  equal  to  zero 
is  surprisingly  good.   It  would  be  useful  therefore  to  be  able  to 
obtain  this  limit  without  solving  the  integral  equation  (9),    for  if 
the  ion  distribution  function  were  more  complicated  this  equation 
might  be  very  difficult.   We  shall  avoid  this  intermediate  step  of 
determining  the  electron  distribution  function  by  showing  that  if 
the  distribution  function  is  factorable  the  electron  pressure  is 
directly  proportional  to  the  number  density  in  the  limit  of  zero 
mass.   Since  the  density  is  a  known  function  of  a  we  may  proceed 
immediately  to  solve  the  field  equation  (l). 

To  show  this  we  return  to  the  definitions  of  the  pressure  and 
number  density 
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P   (a)  =  4  r  f(E,P,Q)  x/ar^  -  (P-a)^  -  q2  dE  dP  dQ 
W(a)   =  -|   r  f  (E,P,Q)  /  VsmE  -  (P-a)2  -  q2  dE  dP  dQ 

and  consider  the  limit  as  m  approaches  zero.   The  range  of  integration 
is  over  those  values  of  E^P  and  Q  for  which  the  square  root  is  real. 
We  observe  that  as  m  approaches  zero  the  range  of  values  that  may  be 
taken  on  by  P  and  Q  approaches  zero.   If  we  assume  that  the  distribution 

function  may  be  factored  so  that  f(E,P,Q)  =  f-j_(P)f2(E,Q)  and  if  f^  is 

12 
a  continuous   function  of  P  then  we  may  evaluate  f^  inside  the 

integrals  at  P  =  a  since  in  this  limited  range  of  integration  the 

difference  between  P  and  a  is  proportional  to  x/iii  .   We  then  have  for 


the  ratio  of  pressure  and  density 

0       (a)  -  /  fo 

•XX         m  J   2 


p    (a)      -   /  f^(E,Q)  (2mE  -  Q^)  dE  dQ 


N(a) 


/  fgCE^Q)  dE  dQ 


where  the  right  hand  side  is  a  constant  with  the  dimensions  of  energy 
and  hence  proportional  to  the  average  energy  per  particle  which  we 
shall  assume  remains  finite  as  m  approaches  zero.   It  is  in  fact  the 
average  energy  per  particle  associated  with  motion  perpendicular  to 
the  magnetic  field. 

If  we  take  f„  ~  5(E-E  )  we  obtain  the  results  of  the  previous 
section  (see  equations  (9)  and  (18)),  namely 


In  the  previous  example  f-,  (P)  is  not  continuous,  but  approaches  a 
continuous  function  as  m  approaches  zero. 
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or 


p    =  2/3  E  N 


P    =  2/3  :^  a 
XX        e 


If  the  plasma  were  Maxwellian,  f _  ~  e  '   ,  we  find 


p    =  W  k  T 

XX 

which  is  the  equation  of  state  of  an  ideal  gas.  As  a  last  example 

let  fp  ~  6(E-E  )  5(Q)  so  that  the  particles  are  confined  to  move 

in  a  plane  perpendicular  to  the  magnetic  field.  The  equation  of 
state  becomes 

p    =  E  N   • 
XX       o 

13 
This  last  result  was  obtained  previously  by  Tonks   although  in  a 

somewhat  different  form.   He  analyzes  in  detail  the  motion  of  the 
particles  in  the  guiding  center  approximation.   It  is  not  surprising 
that  these  two  methods  should  give  the  same  results,  for  letting  the 
mass  approach  zero  is  equivalent  to  letting  the  Larmor  radius  approach 
zero,  or  in  effect  assuming  that  the  Larmor  radius  is  small  in  com- 
parison to  the  distance  within  which  the  magnetic  field  changes  by 
a  significant  fraction  of  itself.   This  is  just  the  basis  of  the 
guiding  center  theory. 

There  is  nothing  in  the  preceeding  analysis  that  restricts  the 
results  to  either  ions  or  electrons.   In  the  transition  layer  problem, 


-^  L.  Tonks,  Phys.  Rev.  97,  1^^1-3  (l955).   See  also  L.  Spitzer,  Jr., 
Astrophys.  J.  II6,  299  (l952). 
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however, the  entire  magnetic  field  variation  takes  place  within  one 
ion  Larmor  diameter  and  so  the  guiding  center  approximation  is  not 
suitable.   The  electron  on  the  other  hand  has  ^3  diameters  in  the 
layer. 
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